p-spin model in finite dimensions and its relation to structural glasses 
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The p-spin spin-glass model has been studied extensively at mean-field level because of the insights 
which it provides into the mode-coupling approach to structural glasses and the nature of the glass 
transition. We demonstrate explicitly that the finite dimensional version of the p-spin model is in 
the same universality class as an Ising spin glass in a magnetic field. Assuming that the droplet 
picture of Ising spin glasses is valid we discuss how this universality may provide insights into why 
structural glasses are either "fragile" or "strong". 

PACS numbers: PACS numbers: 05.50.+q, 75.50Lk, 64.60Cn 



In recent years structural glasses have been studied 
by many researchers using ideas from spin glasses. This 
activity was initiated by Kirkpatrick, Thirumalai and 
Wolynes |l| who observed the similarity between the dy- 
namical behavior of p-spin spin glass models at mean- 
field level and the mode-coupling equation for glasses . 
At mean-field level, the p-spin model has two transitions 
1^ of interest for the structural glass analogy. There is 
a dynamical transition temperature below which the 
spin-spin correlation functions do not decay to zero in the 
long time limit which implies a breakdown of crgodicity. 
It is this transition which corresponds to the transition in 
the mode-coupling equations. The transition at Tc < Td 
is a transition associated with thermodynamic singulari- 
ties (there are none at the dynamical transition Td) and 
below it replica symmetry is broken at the "one-step" 
level It is this transition which is usually identified 
with the glass transition (if any) of structural glasses and 
Tc associated with, say, the Kauzmann temperature Tk 

i- 

It has long been recognised that the dynamical tran- 
sition at Td is an artifact of the mean- field limit ||,§. 
The transition is due to the exponentially large number of 
states which trap the system for exponentially long times 
forbidding the system in the thermodynamic limit from 
reaching equilibrium. However, in finite dimensional sys- 
tems activation processes over finite free-energy barriers 
or nucleation processes will take place and restore ergod- 
icity, a feature which was not incorporated either into 
the early versions of mode-coupling theory. Thus the 
type of transition which occurs at is the only possible 
transition which might exist in finite-dimensional p-spin 
models. In fact in this paper we will argue that this tran- 
sition does not survive either in finite dimensional models 
of p-spin glasses. If the connection between real glasses 
and p-spin models goes beyond the mean-field limit then 
one would expect that there will be no genuine transition 
in real glasses either. 

Our approach will be to start from a commonly used fi- 



nite dimensional version of the p-spin model and study 
it using field-theoretic methods. We find that the field 
theory of the p-spin model (at least when p is three, the 
only value that we shall consider here) is in the same uni- 
versality class as the Ising spin-glass in a magnetic field. 
The possibility of such a connection was suggested by 
Parisi et al. based on their numerical simulations and 
by Drossel et al. |^ from their Migdal-Kadanoff renor- 
malisation group studies, but it has never to our knowl- 
edge been explicitly demonstrated. While both models 
lack time reversal invariance, at mean-field level the p- 
spin model and the Ising spin glass in a magnetic field 
behave quite differently so that their similarity for finite 
dimensional systems cannot be taken as obvious. This 
is because the phase below Tc in the mean-field p-spin 
model has replica symmetry broken at the one-step level 
Q, whereas the transition in the mean-field Ising spin- 
glass in a magnetic field - the transition discovered by 
de Almeida and Thouless (AT) 0, ||] - has the replica 
symmetry broken at all steps. 

Away from the mean-field limit the existence of a tran- 
sition for the Ising spin glass in a magnetic field remains 
controversial. According to the droplet picture , 

the presence of a field removes the phase transition just as 
a field applied to a ferromagnet does. However, advocates 
of the replica symmetry breaking picture of spin glasses 
claim that the AT transition exists in finite dimensions 
|]l2| . In what follows we shall interpret our central result 
of the equivalence of the p-spin model with the Ising spin 
glass in a field according to the droplet model, and not 
discuss the consequences of this equivalence in the replica 
symmetry breaking picture. To support this interpreta- 
tion we have determined the correlation length ^ of the 
p-spin model as a function of the temperature T using 
the Migdal-Kadanoff approximation (MKA) |^ . The cor- 
relation length rises rapidly below a certain temperature 
before saturating at a large value, which is around 30 
lattice spacings in three dimensions, but there is no evi- 
dence within this approach of an actual phase transition 
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and a diverging correlation length. 

In the three-spin model used for the present study, each 
site is occupied by two Ising spins, Oi and r^, and the 
Hamiltonian is 

n{{a},{T})= - ^ (^JV.T.a, + 4'V,T,T, 



<ij> 



7(3) 



(1) 



where < ij > are nearest-neighbor pairs, and the cou- 
plings J^^"^ are chosen independently from a Gaussian 
distribution with zero mean and width J. 

When the signs of all the spins are reversed, the sign 
of the Hamiltonian changes, indicating the violation of 
time reversal symmetry. The thermal averages (ai) and 
(Ti) are non-zero at all finite temperatures. 

In order to derive a field theory for the three-spin 
model we use the usual replica identity for the free energy 
F 



F = -kT\n{Z) 



-kT Van ■ 

n— »0 



1 



(2) 



(n) 

to perform the average over the couplings J^^ , indicated 
by the overline. To evaluate we set up n replicas of 
the partition function Z 



for which a = (3 are to be omitted. We next follow the 
procedure of Pytte and Rudnick |l3[ and define 



a/3 
ij 

3 ^ij 



(6) 



We are left with the following types of sum at each site: 



Tr exp 



Qa/3(i)crfcrf 



V2 



(7) 
tak- 
Re- 



These can be done by expanding the exponential, 
ing the trace and then re-exponentiating the result 
peating this procedure for each site and finally taking the 
continuum limit one obtains, correct to 0((5^^, ^q^), the 
following free-energy density 



+ PcfiPU) - a(,Ql^ + Q*^p + Pip + P:^) 

1 

6 



W2 {PlpQaP + PafiQap) 



(8) 



- Tr n exp [^H {r"}) /kT] . 



(3) 



Tr denotes taking the sum over the values ±1 of the Ising 
spins erf and rf at each site. After bond averaging one 
gets 



Z" = Tr exp 



1 
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^3 '^3 



<i3> 

I ^a. 



, (4) 



where z is the number of neighbors of a site (six for a 
simple cubic lattice) and N is the number of lattice sites. 

The terms involving different sites are decoupled by 
introducing fields q°/^^ and p^^ for each bond < ij > 
using: 



Z" = Tr exp {P^J^nzN) Y[ 



<t3> 



72 a/3 



j2 a/3 
'^P^3 



X exp 

a/3 



/ a/3 a/3=i= . a/3 a/3+ 



./3, 



,/3 



ctB^ a a 
■P^3 ''3^3^3^3 



Prj (fT, < 



3 J 



,(5) 



where all fields have the symmetry Qap = Qpa, ol ^ (3 
and the implied summation convention is used. Notice 
that the free energy derived from the functional T is 
real as can be shown by taking the complex conjugate 
of Eq.(||) and making a variable change P <^ Q in the 
fields, which are of course integrated over to get the free 
energy F . 

The above field theory, while complicated, is of stan- 
dard type. The field theory for the p-spin model obtained 
by Campellone et al. [Q is very unusual and does not 
lend itself to a diagrammatic expansion. Truncating the 
expansion at the cubic terms is adequate provided that 
the phase transition, if any, is continuous. At mean-field 
level, the transition is continuous provided W2 /wi is suf- 
ficiently small (see below). The numerical studies of the 
three-spin model in finite dimensions suggest that the 
transition, if it exists, is continuous, so that we shall fo- 
cus on situations in which the truncation to a cubic field 
theory is sufficient. When this is the case further sim- 
plifications are possible by simplifying the field theory to 
the soft modes associated with the possible transition. 

First of all, it is easy to see from the quadratic terms 
that as the coefficient a is positive the terms involving 
the imaginary parts of Q and P are associated with hard 
modes and so can all be dropped near the putative tran- 
sition. The associated free energy density functional is 
then 



where in the applied summation convention used here 
and in what follows for terms containing a and /3 those 



J" 
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(VA„;3) 



1 



(VP„^) 
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- \w2{AlpB^p + A^fiBlp), (9) 

where A^p is the real part of Qap and B^p is the real part 
of Pap] -t/4: = {r - 2a), wi = {1 + l/\/2)wi, W2 = W2- 
This functional is manifestly real. 

A second simplification comes from observing that 
even at high temperatures {Aap) = {Bap) ^ 0, that is, 
the average of the fields is always non-zero. This is due 
to the term {A'^pBap + AapB'^^). Any transition which 
occurs will involve a breaking of this replica symmetric 
solution. It is convenient to first introduce the new fields 

= {Aap + Bap)/V2 
Dap = {Aap-Bap)/V2. (10) 

Then in terms of these new fields the free energy func- 
tional becomes 

^ = \{^Cap)^ + \{VDap)^ - \t{Clp + Dip) 

- -^'V2{Cap - CapD^p) - -Vi{Cal3Cp~^C^fa (H) 
+ DapDpjCja + DapCpjDja + CapDpjD^a), 

where vi = wi/ \/2 and f 2 = W2 / \/2. If one introduces 
a field Sap such that Cap — Sa/3 + C with < Sap >= 
then one can see that the terms involving the Dap fields 
are hard modes and can therefore be dropped. The free 
energy density functional is then 

^ = -^i'^Cap)'^ - -^tC'lfl ~ -^ViCaljCp^C^a " -V2Cap- 

(12) 

Ferrero et al. jisj wrote down this functional as the 
generic functional for discussing replica symmetry break- 
ing at the one step level. They showed that the transition 
at mean- field level is continuous when V2/V1 < 1. 

The functional in terms of the Sap fields (leaving out 
terms independent of Sap) is 

r = \{VSap)^ -\{t + ViC^ +V2C^)Sa0 

- ^{t + 2v2C)slp-'^ViCSapSaj 

1 1 3 , , 

- QViSaffSpySja " gW2SQ0, (13) 

where in the term ■^viCsapSa-y, P ^ 1- Writing the 
quadratic terms in the form 

-^(ySafi)'^ + Map,'i5Sa0SyS, (14) 

it turns out that the matrix M has three distinct eigen- 
values . The limit of stability of the replica symmetric 
high-temperature phase at mean-field level is determined 



by the vanishing of the smallest of these. AT showed 
that the corresponding eigenvectors span an ^n(n — 3) 
dimensional subspace of the ^nin — 1) dimensional space 
defined by the Sap- This is called the "replicon subspace" 
p6[ and since only the replicon modes go soft at the AT 
or replica symmetry breaking transition it suffices again 
to retain only these modes for a study of the critical be- 
havior of this transition. It is useful to introduce the 
projection operator Sap.-iS which projects any field Sap 
onto the replicon subspace Q. Then qap = Saf3,jSS-yS, 
where the sum is over distinct pairs (7(5). The matrix 
elements are Sap.ap — SjSap.a-y — 3/2, (/3 ^ 7) and 
Saff.-yS = Ij (0)/3)7)<5 all different) when n is set to zero. 
Then our final free energy density functional is 

\C^qo^fi)^ + ^iQlp-^Viqaf3qi3jqya-^V2qli3^ (15) 

where the coefficient i is a measure of the distance from 
the transition (if any). Bray and Roberts 1^ first ob- 
tained this functional in their study of the critical expo- 
nents in dimension 6 — e at the AT line. 

Thus the ultimate field theory for the transition (if 
any) in the p-spin model is that of the Ising spin glass 
in a magnetic field. Bray and Roberts |^ obtained the 
recursion equations for the coupling constants vi and V2 
correct to order e but did not find any stable fixed point. 
They suggested that a possible explanation for this might 
be that there is no transition. We next provide some ev- 
idence using the Migdal-Kadanoff renormalisation group 
procedure that this indeed may be the correct interpre- 
tation of their perturbative renormalisation group calcu- 
lation. 

The MKA is a real-space renormalization method 
which we have previously applied |^ to both the three- 
spin model and the Ising spin glass in a field. For both 

(71) 

models the behavior of the couplings such as J^^ un- 
der the renormalisation group was to fiow to the high- 
temperature sink i.e. zero, implying that there was no 
transition. We have now investigated how the effective 
couplings decrease as a function of the lengthscale L and 
find that the variance of Jj'"' decreases as exp(— L/^), 
where ^(T) is the correlation length of the system. The 
dependence of 5(T) on temperature is shown in Figure |l| 
and it is obvious that because of the rapid rise of ^(T) 
when T « 2 J Monte Carlo numerical studies on systems 
of linear dimensions less than 10 could easily appear to 
have a phase transition |^ . 

What can we deduce from these results on p-spin 
glasses for the physics of structural glasses? The connec- 
tion of the transition at Tc in the p-spin model (if any) 
with the structural glass transition at around Tk (if any) 
has been much discussed |^ . The relationship can really 
only go beyond an analogy if the two systems belong to 
the same universality class like, for example, the Ising fer- 
romagnet and the liquid-gas critical point transition. At 
first sight any relationship between the p-spin model and 
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T/J 

FIG. 1: The correlation length of the three-spin model in 
three dimensions as function of T/J, as obtained using the 
MKA. 



is not large but instead of order one. It is our contention 
that the structural glass problem is similar to the Ising 
spin glass in a field, and if that field is large the correla- 
tion length is small |^ and the non-Arrhenius behavior 
produced by a rapidly increasing correlation length will 
be absent in those circumstances. The validity of this 
picture would be greatly strengthened if experimental 
evidence of a growing correlation length at temperatures 
around Tk in fragile glasses were to be experimentally 
observed. 

One of us (M.A.M) would like to thank A. J. Bray 
and A. Cavagna for numerous discussions. B.D. was sup- 
ported by the Deutsche Forschungsgemeinschaft (DFG) 
under Contract No Dr300-2/l. 



phenomena in supercooled liquids seems unlikely. The 
p-spin model has quenched disorder, and there is no dis- 
order at all in the glass transition problem. If there is 
a connection it would probably arise from similarities of 
the pure states in the two systems. The one-step replica 
symmetry broken state below Tc implies that there are 
many pure states present each with the same overlap with 
each other |^ . If a structural glass is formed via a phase 
transition it might have pure states with similar proper- 
ties and then the functional of Eq. (|l^) might describe 
its transition correctly. 

It is our contention though that there is no phase tran- 
sition in the finite-dimensional p-spin model (in any finite 
dimension) and hence discussions of universality classes 
are not possible. However, the fact that a large corre- 
lation length can arise at low temperatures means that 
scaling and universality ideas should have a degree of util- 
ity and hence that the p-spin model may provide a qual- 
itative model for structural glasses. According to this 
possibility one could argue that the more rapid growth 
than Arrhenius for the relaxation time r seen in some 
glasses ( the "fragile" glasses) is because the relaxation 
time grows as 

T = ToCxp[Be(T)'^/fcr] (16) 

accepting the connection to the droplet model of spin 
glasses, where tq is a microscopic relaxation time, and 
the constant B would be expected to be of order J. The 
exponent ^ would be that of the zero field Ising spin glass 
pO| . The rapid rise of S,{T) depicted in Figure [I] would 
then produce a very non-Arrhenius temperature depen- 
dence of the relaxation time r. On the other hand some 
glasses (the "strong" glasses) do have an Arrhenius tem- 
perature dependence of their relaxation time. We would 
interpret this as arising when ^(T) at low temperatures 
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